Coupling bosonic modes with a qubit: entanglement dynamics at zero and finite-temperatures 
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We consider a system of two iso-spectral bosonic modes coupled with a single two-level systems i.e., a qubit. 
The dynamics is described by a mode-symmetric two-modes Jaynes-Cummings. The entanglement, induced 
between the two bosonic modes, is analyzed and quantified by negativity. We computed the time evolution 
of negativity starting from an initial thermal state of the bosonic sector for both zero and finite temperature. 
We also studied the entangling power of the interaction as a function of mode-qubit detuning and its resilience 
against temperature increase. Finally a two-qubit gates based on bosonic virtual subsystem is discussed. 



INTRODUCTION 

In several of the existing implementation proposals for 
quantum information processing 1 1 ] hybrid system are in- 
volved. By hybrid here we simply mean that both discrete i.e., 
atomic, and continuous e.g., vibrational modes in an ion-trap 
0, degrees of freedom are present and interacting. Typically 
the former are chosen to play the role of quantum informa- 
tion carrier whereas the latter are treated either as informa- 
tion buses (to couple on demand remote discrete systems e.g., 
cavity modes coupling atoms 01) or are hold responsible 
for decoherence effects e.g., phonons in quantum dots. In ei- 
ther cases the attention is, in a sense, mostly focused on the 
information-carrying discrete sector, and the bosonic modes 
are eventually traced away. 

In this paper we would like, so to speak, to reverse this logic 
and to address the problem of bosonic modes interacting using 
discrete systems as interaction medium. We will consider the 
simplest instance of this situation: two iso-spectral bosonic 
modes coupled through the interplay of one qubit. We con- 
sider here the bosonic degrees of freedom as the relevant ones 
for encoding information, while the qubit degrees of freedom 
will be traced out. Information will be encoded in a finite- 
dimensional subspace of the two-mode bosonic Fock space. 

Of course the use of continuous variables i.e., light modes 
in the context of quantum information processing has been al- 
ready massively investigated in the recent literature. On one 
hand there is the effort to handle with the infinite dimensional 
Fock space of bosonic degrees of freedom by mapping it onto 
finite dimensional ones by choosing ad hoc states. On the 
other hand a careful choice of a suitable subspace of the Fock 
space can also provide a proper way for encoding and process- 
ing qubits 1 5] The first approach is the one contained e.g. in 
0] in which the gaussian states are involved in view of their 
key feature of being easily mappable onto the finite dimen- 
sional space covariance matrices. The second approach in- 
stead aims at finding information-encoding states whose time 
evolutions remains inside a finite dimensional subspace, at 
least for properly selected time instants 0,0,0]. 

This second approach is the one followed in this paper. 
Given the interaction Hamiltonian, a proper encoding i.e., a 
two-dimensional Fock subspace, is searched. The require- 



ment that one can individuate time instants in which a non- 
trivial unitary transformation of the encoding subspace is en- 
acted. 

More specifically, in this paper we shall analyze a sym- 
metric two-modes Jaynes-Cummings model. The entangle- 
ment properties of the bosonic subsystem are investigated by 
computing the dynamical evolution of the negativity 1 10] of a 
generic state of the system. Analytical and numerical results 
respectively are found for the vacuum and the thermal initial 
state. The latter case is also studied to prove the entanglement 
resilience against temperature and qubit-modes detuning. The 
entanglement capabilities of our system are then analyzed by 
introducing a convenient notion of entangling power II lUl 211 
and by observing its behaviour as a function of the model pa- 
rameters. Finally a unitary quantum gate for this kind of inter- 
action, is contrived. The realization of this gate however will 
require the use of virtual subsystems 1 13] and of the related 
entanglement relativity concept 1 14]. 



THE MODEL 

The physical system studied in this paper is the one de- 
scribed by the two-mode Jaynes-Cummings hamiltonian |7, 
0: 
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H = cr z e + y] hujiblbi + yVo-_&j + a+bi) (1) 
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The symbols bi and b\ stand for annihilation and creation 
of a particle (photon or phonon) in the mode i, while 7 is 
the coupling constant and e is the energy difference between 
the two states. The hamiltonian simplifies when i) the two 
modes have the same energy uji = 0J2 = u), and ii) the system 
achieves resonance e = fku. For the sake of simplicity we deal 
first with this case. In order to diagonalize the hamiltonian we 
perform a unitary transformation of the two mode operators, 
introducing two new modes: b + = -^(bi + 62) and = 

^(61—62), as well as their conjugated 6^, and the inverse 



equations: b\ = (b + + 6_) and b 



75^ 



b-). The 
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hamiltonian then becomes: 

H = cr z e + frhjtf + b + + j(cr-b\_ + cr+b + ) + twb[_b- (2) 

As we can see, now the hamiltonian is the single-mode 
Jaynes-Cummings hamiltonian plus a commuting term (the 
last one). The interaction involves now only one mode while 
the other does not evolve. The evolution of a genericjpure) 
state is now easy to obtain following for example ref. 1l7ill 811 : 
if the system is prepared in the initial state with the qubit 
in its excited state (e) and the two bosonic modes occupied 
respectively with m and n 2 particles (|^) = \e, 77,1,77,2) = 
E n+ , tt _ a„ +! „_ |e, n+, ti_)), the evolution will be: 

= ^2[A ni ,n 2 {t)\e,ni,n 2 ) 

ni,n 2 

+ B nit7l2 (t)\g,ni + l,n 2 ) + C ni , n2 (t)\g,ni,n 2 + 1)] 

n_)_ ,n_ 

- i,s„ + (i)| i g,n + + l,n_)] (3) 



ENTANGLEMENT DYNAMICS 

In this section we will study the bosonic entanglement dy- 
namics enacted by Q; we will start from different initial 
preparations i.e., zero and finite temperature. As entangle- 
ment measure we adopt negativity introduced in Ref. ficll 

Vacuum State 

If the system is initially in the ground state with respect 
to the two bosonic modes, it is straightforward to obtain: 

Pm,m = Pn+,n- = |e, 0, 0) (e, 0, 0| being = n± = 0. 
The time evolution of this state is given by: 

p(t) = ^(t)|e ) 0,0)±(e ) 0,0|+ fl g(t)| fl ,l,0)±< fl ,l,0| 
+ ico(t)so(t)(|e,0,0) ± ( 5) l,0| + | 5 ,l,0)±(e,0,0|) 

(6) 



where c n (t) = cos £l n t, s n (t) = smfl n t and 0„ = js/n + 1 
is the (half) Rabi frequency. The coefficients a n+ ,n_ are ob- 
tained by writing the initial state as a sum of states with fixed 
number of particles in the new modes: 

1 t ni t n2 

|e,m,n 2 ) = , , , b\ b 2 |e,0,0) 



\fni\n 2 \ 

^ \ rii+n 2 ^ 
ni+Tt2 



Here and in the following we denote with |)i 2 and \)± the 
elements of the two basis sets {|?ii)8)|n 2 )} and 
respectively. 

We choose to use the negativity as a measure of entan- 
glement for bipartite states. The negativity is defined as the 
sum of the negative eigenvalues of the partial transpose (trans- 
posed with respect of one of the two subsystems) [19, 20] of 
the density matrix of the state. In order to compute the nega- 
(b\_ + b^ ) ni (&t — &t )™ 2 |e, 0, 0) tivity of the state with respect to the bipartition defined by the 

two bosonic modes we trace out the qubit degrees of freedom: 
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e,0,0) 



n+ ,n_ 



a 



n_|_,n_ |c, Tl^- j Tl—J 



(4) 



having identified, in the last passage, i+j and 77,1 + 77 2 — i — j 
with 77 + and 77_ respectively. 

With this formalism we can in principle describe the evolu- 
tion of a generic density operator: 



pit) = cg(i)|0,0)±(0,0| + «g(t)6! h |0,0) ± (0,0|6 + 

+ c 2( t )| 0) 0) 12 (0,0| + ^(6l + 4)|0,0)i 2 (0 ) 0|(6 1 + & 2 ) 

= cg(t)|0,0) 12 <0,0| + ^(|l,0) 12 <l,0| 
+ |1,0) 12 <0,1| + |0,1) 12 (1,0| + |0,1) 12 <0,1|) (7) 

The reduced density operator can be now represented as a 
matrix in the |m) (8 \n 2 ) basis, and its partial transpose p Tl 
with respect to (e.g.) the subsystem of mode 1 reads: 



P / y Pn 1 ,n 2 Pn 1 ,ri2 



«1 ,"2 

E 



|e,ni,n 2 )(e,77,i,n 2 | 



til ,n.2 



Vn\,n 2 X 
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(8) 
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^2 a «+,™- a <,n'„ l e ' n +' n -)( e : n '+> n '- 



The negative part of the spectrum of this matrix gives the 
analytical result for the negativity: 



(5) 

In general this straightforward procedure can be quite difficult 
to handle because of the big amount of terms involved in the 
sum. However in two cases the expression of the time evolu- 
tion of the initial state highly simplifies. They are the vacuum 
and the thermal state cases. 



(9) 



whose graphical plot is displayed in fig.Q. It is interesting to 
note the periodic behaviour of the negativity which vanishes 
with period 77/7. On the other hand we can see that this kind 
of interaction is able to entangle very much the system, giving 
raise to maximally entangled states with the same period 7T/7. 
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In fact as t = tt/2j + K(ir/j) the state vector is (except for 
a phase factor): 



-L(|i,o) 12 + |o,i) 12 ) 



(10) 



In fig.Q, besides the negativity, we plotted the linearized Von 
Neumann entropy too: Sl(p) = 1 — 2>i(T>2p) 2 , with Tri 
meaning trace over the subsystem i. Its analytical expressions 
for the present state is given by: 



Si = 1 - ck - 



(ID 



As long as p(t) is a pure state, the Von Neumann entropy is a 
measure of entanglement and its maximum coincides with the 
negativity of the maximally entangled state. The state p(t) is 
pure when t = ir/2-f + K(n/2'y), as it is easy to recognize 
by computing its spectrum: er[p(<)] = [0, 0, Cg(t), Sg(t)]. It is 
worth noting that when the state is no more pure the Von Neu- 
mann entropy may exceed the maximum entanglement value 
for pure states, then being a no more reliable measure of en- 
tanglement, as fig.© shows. 




yt 



FIG. 1 : Negativity as a function of time for the vacuum state (solid 
line). The linearized Von Neumann entropy is plotted in dots for 
comparison: the two lines touch each other when p is a pure state. 

This result holds even for the more general case of a multi- 
mode hamiltonian, and provides a quick way to generate the 
so-called W n-qubit state 

|^) = -L(|l,..,0)+... + |0,..,l)) (12) 
Jn 



This feature will be exploited in the design of a quantum gate 
in the following. 



Thermal State 

Let us move now to the case in which both bosonic modes 
are in a thermal state. The bosonic system will be described at 
the initial time as a tensor product of the two thermal states: 
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-j3hujni 



-f3huj(ni-\-ri2) 



-{3Huj(n+-\-Ti- 



(13) 

where, as usual, (3 = -^t an< ^ n * s usec ' m tne °P~ 
eratorial meaning (n = b'b). The last identity of eq. 
Jl 3I > allows us to write the density operator directly in the 

|n+) ® |n_) basis: p = £ n p ni .n 2 K, n 2 )(n 1 , n 2 \ = 

2n+,n_ Pn + ,n-\ n +i n -){ n +i n -\ where the sum 8 oes as 
usual from zero to infinity. The coefficients p n account for the 

Bose-Einstein probability distribution, according to which: 



Pn 



and the mean boson number is 

(n) = — 



(14) 



(15) 



The time evolution of this thermal state, after tracing out the 
qubit, will be: 

p(t) = Pn+,n-[($ l+ {t)\n+,n-){n + ,n-\ + 

n+ ,71_ 

+ 4 + (i)|n + + l,n_)(n+ + l,n_|] (16) 

The matrix form of this operator, of course, is an infinite 
square matrix whose spectrum can be computed numerically 
thanks to the sparse structure of the matrix itself. By increas- 
ing the temperature the higher energy levels begin to con- 
tribute to the density and the infinite sum of terms has to be 
truncated carefully, taking into account more and more levels. 
The matrix elements of p Tl can be written as a sum of matri- 
ces with fixed number of particles weighted with their thermal 
probability distribution. By combining eq.© and © as well 
as their inverse, it is possible, after some tedious but straight- 
forward calculations, to write down the matrix elements of 
p Tl as: 
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# = A 



S d' x -1 (*) ^ y^^^^^Y^^TC^y ^ _ X ) fe+J Pd' 1 - lPdi ^,n+2+( AT+2) (i+j) - (fe+O ^,1^+2+ ( JV+2) (fc+i ) - (i+J ) ] ( 17 > 



where d 2 = n — d\\ d'% = n + 1 — d[ and i,j = 0, d\\ 
k,l = 0, d 2 . 

The matrix has a multi-diagonal form, meaning that the 
non-zero elements are placed in diagonal lines parallel to 
principal diagonal of the matrix. The computation of its 
spectrum and in particular of its negative part (negativity) 
has been done in a numerical way. The series appearing in 
ea.in\ has been truncated using as a convergence criterion: 
Sn=o J2d 1= oPdiPd 2 ~ !■ The results are shown in fig.0 in 
which the time evolution of the negativity is plotted for dif- 
ferent values of the energy/temperature ratio r\ = -^Hp . As 
one can expect, the negativity peaks get lower as tempera- 
ture raises. Nevertheless, surprisingly enough, the negativity 
remains different from zero for long time even at high temper- 
ature. At T — we recover the same results of the previous 
subsection. 
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FIG. 2: The negativity behaviour of the thermal state as a function of 
time is plotted for four different temperature. The latter are labelled 
with different values of the parameter I/77: for the solid thin gray 
line, 0.5 for the dashed line, 1 for the dotted line, 5 for the solid black 
line. 

In order to make our analysis slightly more general, 
we now extend these results to the case in which the 
qubit and the bosons are off -resonance. This occurs when 
the qubit energy splitting and the bosonic mode frequen- 
cies are different. In this situation the coefficients c n 
and s n have a more complicated time dependence fnll : 



c„ = [cos(fi n i) — sin(fi n t)] exp(iAt/2) and s n — 
(ijy/n + l/fln) sin(fi n i) exp(z A£/2), while the (half) Rabi 
frequency becomes: il n = ^y/ A 2 + A^f 2 {n + 1), with A = 
(e — huj)/h. Fig.(|3} shows the negativity behaviour at T = 
for increasing values of A, while fig.© shows how the neg- 
ativity varies in time for different temperatures with a fixed 
non-vanishing value of the detuning (A = 1). As a general 
remark, it is worth noting that an increasing detuning causes 
the peaks of the negativity to lower at fixed temperature. In a 
similar way the high temperatures lower the negativity as al- 
ready known from the resonance case (fig.(|3 and @ are quite 
similar). Finally for high values of A the negativity oscillates 
with frequency 7 2 /A as a result of a second order effective 
interaction hamiltonian H e ff\ 

H eff = m ef fb\b 2 + h.c. (18) 

where fi e // = 7 2 /A. In fig.O) the dotted line plots the neg- 
ativity for A — 2 and 7 = 1. Since A >> 7 we can con- 
sider H e f f as a perturbation and we can estimate the ener- 
gies of the evolving states as the unperturbed energy splitting 
e ~ A plus a second order energy correction 7 2 /A. The 
same result is obtained by considering the original hamil- 
tonian with 2«o = y/A 2 + 4 7 2 = Ay/1 + 4 7 2 /A 2 ~ 
A[l + 2(7 2 /A 2 )] = A + 27 2 /A. The negativity oscillates 
with double frequency with respect to the state vectors, since 
its time evolution is dictated by Sq and c 2 . Therefore its ap- 
proximated Rabi frequency will be 2f2rj = A + 2 r y 2 / A = 3 
as it is clearly shown in fig. (0. 

ENTANGLING POWER 

The above qualitative remarks have been quantitatively 
proven by defining the entangling power lfl2ll of the hamil- 
tonian as: 

E P (T)=supAf[p T (t)} (19) 

where r is the time variation period of p, and has been cho- 
sen long enough to allow all the frequencies contributing the 
time evolution to be taken into account. The reference state p 
is a thermal state and the entangling power is computed as a 
function of the temperature. Fig.Q shows its behaviour for a 
range of the ratio 1/ry going from to 10. Although the en- 
tangling power quickly decreases as temperature raises, it is 
however different from zero for a wide range of temperatures. 
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FIG. 3: Time variation of the negativity at T = for different values 
of the detuning parameter: A = (thin gray line), A = 1 (dashed 
line), A = 2 (dotted line). Here the coupling constant is set as 7 = 1. 



FIG. 5: The entangling power of the hamiltonian for the thermal 
state as a function of the temperature. As long as I/77 is lower than 
1, Ep decreases very quickly, then it has a slow variation, and for 
high values of the parameter it remains almost constant. 





FIG. 6: Entangling power over all the off-resonant hamiltonian la- 
beled by the detuning parameter A going from to 10. 



FIG. 4: Plot of the negativity time dependance in the case a fixed 
detuning A = 1. Three different values of 1/rj have been chosen: 
(thin gray line), 0.5 (dashed line), 1 (dotted line), 5 (solid line). As 
in the previous case: 7 = 1. 



The reference state is again the thermal state, and the 3- 
dimensional plot of Ep is represented in fig.Ql. The most 
interesting feature of Ep is its monotonicity with respect to 
both T and A. 



As a second step we studied the entangling power of a wider 
set of hamiltonian. As done in the previous section, we intro- 
duce the detuning parameter A in the hamiltonian and redefine 
the entangling power as: 

Ep(A) = supjV [p A (t)l (20) 

wherep is now the vacuum state, chosen as reference state. 
The result is displayed in fig.© in which we see the entan- 
gling power quickly decreasing down to zero, as A increases. 

In order to combine the two previous results, we finally de- 
fine the entangling power as: 

E P (T, A) = sup N[p${t)] (21) 
ter 



BUILDING A QUANTUM GATE 

In this section we will examine the possibility of using 
the above analyzed entangling capabilities to build a quan- 
tum gate. The idea is to encode quantum information into 
the bosonic degrees of freedom by selecting a suitable finite- 
dimensional subspace of the Fock space. As mentioned in the 
introduction one would like to see whether with the proper en- 
coding one can find out specific operating times such that the 
dynamics enacted by the Hamiltonian {1} amounts to a non- 
trivial transformation of the encoding subspace while acting 
trivially in the qubit factor. This latter requirement being of 
course due to the necessity of avoiding entanglement between 
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FIG. 7: Three-dimensional plot of the entangling power for both A 
and \ jr\ going from to 5. It is clear the monotone character of Ep. 



the bosons and the qubit that would in turn result into deco- 
herence. 

We consider as the encoded qubit the first two states of 
harmonic oscillator of the modes of the system studied so 
far. The logical computational basis is built by identifying: 
|0)i = I <7, 0) 12 and 1 1) ^ = \g, l)i2- This is in a sense the 
most natural basis choice. At variance with the previous cal- 
culations we prepare the system with the qubit in its ground 
state 1 18]. The two qubits evolution can be studied within the 
following scheme: 



|00) L = 


l5,0,0)i 2 


|10>L = 


l5,l,0)i 2 


|01)l = 


M,l>12 


|h)l = 


Is, 1,1)12 



The time dependence of the four basis states is given by: 



|<7,0,0) 12 = | 5> 0,0>±->|5,0,0>± 

Iff,l,0)i2 - -±=(| 5 ,l,0)± + | ff ,0,l) ± ) 

Is, 0,1)12 = -j=(\ g ,i,o)±-\g,o,i)±) 

I.9,l,l>i2 - -j=(\g,2,Q) ± -\g,0,2)±) 



-j=(co(t)\g, 1, 0)± - is (t)\e, 0, 0))± + ±=\g, 0, 1)± 
-j=(ca(t)\g, 1, 0)± - is (t)\e, 0, 0)±) - -j=\g, 0, 1)± 
-j=(cx(t)\g, 2, 0)± - isx{t)\e, 1, 0)±) - -^=\g, 0, 2)± (23) 
I 



With respect to the case in which the qubit was initially in its 
excited state, in this case the global state of the system evolves 
according to c, i+ _i and s n+ -x instead of c n+ and s n+ . 

Given this evolution table, it is possible to design a non- 
trivial unitary transformation on the computational basis by 
properly tuning the frequencies appearing int eq. i23\ . By 
inspecting at the dynamical evolution of the basis it is clear 
that one can enact the following gate: 



U = 



Z 1 







-l 



o 



1/ 



(24) 



Indeed, it is sufficient to find a time instant t g for which 
co(t g ) = -l,ci(t g ) = 1 and s (t g ) = si(t g ) = to real- 
ize the above unitary transformation. These conditions can be 
achieved by choosing the variables t g ,j,A such as the fol- 
lowing equation is fulfilled: 

n _ 2N + 1 

fh ~ 2M 

where N, M are natural numbers, that means that for flo, fii 
we can choose Vgo, 9i such that qo/qi G Q. Thus once we 



(25) 



chose qo,Qi, we get: 

7 
A 

to 



2^ ql 



% 



fx 



■^-(2N + 1) = ^—2M 

2<l0 2?1 



(26) 



whence we must add the conditions qx > qo and 2q$ — qi>0 
in order to have 7, A G K. As an example consider qo = 
Q = 3 and qx = fii = 4; then Aj 2 = 7, A 2 = 2 and 
t g = 7T + 2nK. 

The above sketched quantum gate however is the well 
known swap gate together with a single qubit 7r-phase shift. 
Unfortunately it is well known as well that this kind of gate, 
unless ancillary qubits are introduced, is not an entangling 
one. Anyway it is clear that no other unitary transformation 
is realizable within the proposed scheme. In fact as soon as 
t 7^ the evolution brings all the basis states out of the com- 
putational space, except for the first one. The fourth state 
comes back into the computational space only when ci = 1 
and sx = giving raise to the unchanged initial state. The 
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same goes for the the second and the third (cq = 1 and 
si = 0), which have however the extra possibility of getting 
(-)exchanged (co = — 1 and s = 0). The situation would not 
change even if we wanted to start with the qubit in its excited 
state. Again we would have to require that all c„ = e ±I( ^ and 
s n = 0, where the minus sign is suitable only for the second 
and third state. Thus we could just have the same phase shift 
for all the states. This is due to the very fact that some of 
the coefficients (and so the frequencies) appear in every ba- 
sis state, preventing the latter to dephase one from another. 
So one possibility to build a real quantum gate is to choose 
a different computational basis which, although it may not be 
scalable, allows however each state to evolve according to fre- 
quencies different from the others. 

We can in principle build a phase-shifter gate by mapping 
the physical basis on the logical basis in a not straightfor- 
ward way. The previous exercise and eq. Q revealed that 
in order to have the basis last state's phase shifted alone, 
we need the evolution frequencies of this state to be differ- 
ent from the other states frequencies, at variance with what 
happens e.g. in eq. i23\ . A different choice of the ba- 
sis would allows us to overcome this obstacle. One should 
choose the computational basis checking that: i) the states 
\i>) x = ±(N l \X,a x ,b x )± + ... + N n \X,y x ,z x )±) have 



we can let them evolve according to eqs. d23l > and J26l >. and 

get: 



^ a x ',b x >,...,y x ',z x >, with x,x' = 1, 



and X = e,g. ii) one is able to find £ 5 ,7, A such that 
il ax , f2 Zx are commensurable frequencies. Unfortunately 
one has to deal with several commensurability conditions that 
are not obvious to be satisfied. Therefore we will not further 
consider this procedure any but we will focus on a second pos- 
sibility, which is somehow more elegant. 

We can choose a virtual bipartition fTill such that, with re- 
spect to this new bipartition, the swap gate becomes an en- 
tangling one. It is known indeed that the entanglement prop- 
erty of a state must be related to the choice of the subsystems 
between which the quantum correlation is measured 1 14l l2 ill . 
The swap operator itself can be an entangling one if we change 
the subsystems (the qubits in our case), by redefining the 
modes and the basis states. We know for example that the ac- 
tion of the swap operator on the Bell basis states induces a sign 
change of the antisymmetric state. Then we can change the bi- 
partition in such a way that the new subsystems (qubits) are 
no more the two oscillation modes 12111 . but \x) and |A) where 
|X>®|A> = IX A > = !*(*)*) = ^(|0,0(1))±|1, 1(0))) with 
X = \& and A = +, — . Thus by setting the computational 
basis as the Bell one and identifying: 



|00)i 
\01)l 
|10)l 

|H>L 



1 

71 
1 

71 
i 

71 
i 

71 



(| 5 ,0,0) 12 + | ff ,l,l) 12 ) 
(Is, 0,0) 12 -|5, 1,1)12) 



(\g, 0,l)i2 + |s, 1,0)12) 



U' = 



/ 1 \ 
10 
10 

\ -1 / 



(28) 



(Is, 0,1) 



15,1,0)12) 



(27) 



Please note that, with respect to the previous case, now we 
are not simply changing the basis, we are indeed changing the 
computational subsystems. With this new choice the operator 
(1281 is a controlled 7r-phase shift which can in principle be 
used for quantum computation tasks. 

It is important to stress that the use of virtual bi-partitions 
makes now a non-trivial task the realization of single qubit 
gates i.e., those with a non-trivial action on just one of the 
(virtual) subsystems 12111 . On the other hand the aim of this 
section was to discuss some issues related to gate building in 
the context of the dynamics Q rather than proposing a uni- 
versal set of gates. For this latter more ambitious goal, it is 
well known that some sort of non-linearity, possibly effective 
or even measurement-induced, between the bosonic modes is 
necessary [22]. 

CONCLUSIONS 

In this paper we have studied the behaviour of the negativity 
of the quantum states of two iso-spectral bosonic modes sym- 
metrically coupled with a common two-level system. By trac- 
ing over this latter we showed that the interaction creates max- 
imally entangled states for T = in short time (tt/2'j). We 
then turned to study how entanglement production depends 
on the temperature by considering as initial bosonic state the 
thermal one. Not surprisingly the more the temperature raises 
the more entanglement decreases. Nevertheless, we observed 
how entanglement resists even up to relatively high T. A neg- 
ativity peaks damping is also observed as the result of a pro- 
gressive increase of the detuning of two-level system and the 
frequency of the bosonic modes: the more the system departs 
from the resonance condition i.e., detuning zero, the lower 
maximally achievable negativity gets. 

The analysis of the entangling power of the interaction sug- 
gested us the way of designing a quantum gate adopting as 
computational basis the degrees of freedom of the two bosonic 
subsystems (the modes). The choice of the proper computa- 
tional basis however turned up to be not straightforward, in- 
volving the introduction of virtual subsystems in order to be 
able to achieve a non-vanishing entangling power. We would 
like now to conclude by briefly commenting upon possible 
further investigations suggested by the study reported in this 
paper. 

It is clear that one can consider the possibility of imple- 
menting the theoretical framework discussed in the above 
with systems different from the light-matter interaction usu- 
ally chosen as the Jaynes-Cummings hamiltonian reference 
model. By properly adjusting the parameters A, 7 one might 
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wonder whether, for example, is possible to set up a con- 
trollable electron-phonon interaction in a semiconductor het- 
erostructure. Optical phonons can interact via a Frolich cou- 
pling with electrons in the valence band of a semiconductor 
superlattice designed to have only two sublevels. With re- 
spect to the present model the latter would involve more then 
one qubit i.e. atoms, of course. This would force one to study 
a generalization of the simple model discussed above. For 
instance one might have to consider some tailored periodic 
potential profile instead of a single two-level step. 

The study of electron-mediated quantum correlations be- 
tween phononic modes and their potential usage of quantum- 
information carriers is a particularly good illustration of the 
sort of logic which motivated this paper (see introduction). In- 
deed phonons represent an almost prototypical example of an 
incoherent decoherence-inducing system whereas electronic 
degrees of freedom are typically the ones proposed for the en- 
coding the quantum information fisHlfill . Phonons, typically 
in a thermal state, are traced out while ways to enact coher- 
ent manipulations of electrons are contrived. In the future we 
would like to investigate whether an exchange of these roles 
of fermionic and bosonic degrees of freedom can be usefully 
figured out. 
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